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A  numerical  study  employing  discontinuous  Galerkin  method  demonstrating  net  surface  heat 
reduction  for  a  cylindrical  body  in  Mach  17  hypersonic  flow  is  presented.  This  application 
focuses  on  using  sinusoidal  dielectric  barrier  discharge  plasma  actuators  to  inject  momentum 
near  the  stagnation  point.  A  5  species  flnite  rate  air  chemistry  model  completes  the  picture  by 
analyzing  the  effect  of  the  actuator  on  the  flow  chemistry.  With  low  velocity  near  the  stagnation 
point,  the  plasma  actuator  sufficiently  modifles  the  fluid  momentum.  This  results  in  redistribution 
of  the  integrated  surface  heating  load  on  the  body.  Speciflcally,  a  particular  configuration  of 
normally  pinching  plasma  actuation  is  predicted  to  reduce  the  surface  heat  flux  at  the  stagnation 
point.  An  average  reduction  of  0.246%  for  the  integrated  and  a  maximum  reduction  of  7.68%  are 
reported  for  the  surface  heat  flux.  The  temperature  contours  in  the  fluid  flow  (with  maximum 
temperature  over  12  000  K)  are  pinched  away  from  the  stagnation  point,  thus  resulting  in  reduced 
thermal  load.  Plasma  actuation  in  this  configuration  also  affects  the  species  concentration 
distribution  near  the  wall,  in  addition  to  the  temperature  gradient.  The  combined  effect  of  both, 
thus  results  in  an  average  reduction  of  0.0986%  and  a  maximum  reduction  of  4.04%  for  non¬ 
equilibrium  calculations.  Thus,  this  study  successfully  demonstrates  the  impact  of  sinusoidal 
dielectric  barrier  discharge  plasma  actuation  on  the  reduction  of  thermal  load  on  a  hypersonic 
body.  ©2014  AIP  Publishing  LLC.  [http://dx.doi.Org/10.1063/l.4898862] 


1.  INTRODUCTION 

Active  flow  control  is  a  topic  of  great  interest  due  to  its 
vast  practical  applications.  Amongst  various  flow  control 
methods,  plasma  devices  have  many  positive  characteris¬ 
tics.^  They  have  been  applied  from  subsonic  to  supersonic  to 
hypersonic  flow  regimes.  Specifically,  extensive  studies  have 
focused  on  plasma  control  of  subsonic  flows. A  good 
amount  of  work  has  also  been  done  in  the  supersonic  and 
hypersonic  regimes^”^^  with  the  aim  of  reducing  net  mechan¬ 
ical  and  thermal  loads  on  the  vehicle.  However,  flow  control 
methods  at  higher  speeds  have  significant  energy  budget 
penalty.  As  an  alternative,  we  propose  the  use  of  micro¬ 
second  pulsed  Dielectric  Barrier  Discharge  (DBD)  plasma 
actuator^’^^^^  for  the  net  reduction  of  thermal  load  on  a 
space  vehicle  in  hypersonic  re-entry,  with  relatively  lower 
amount  of  power  requirement. 

Plasma  actuators  have  many  advantages  including 
design  simplicity,  rapid  response,  and  lack  of  moving  com¬ 
ponents,  surface  compliance,  and  possibility  of  application 
to  the  most  receptive  location,  relatively  low  power  budget, 
and  robustness. They  can  be  categorized  into  non-thermal 
and  thermal  plasmas.  Non-thermal  plasmas  are  known  to 
induce  reasonable  air  flow  speeds  of  1  10  m/s  within  a 
millimeter  height  from  the  wall  in  quiescent  flow  condi¬ 
tions.  Thus,  they  provide  effective  control  for  free  stream 
velocities  ranging  from  low  subsonic  to  supersonic 
flows. Additionally,  they  have  low  energy  costs,  since 
majority  of  the  electric  energy  goes  into  the  kinetic  energy 


of  the  electrons  rather  than  heating  of  the  surrounding 
gas.^"^  The  examples  of  these  include  corona  discharge  and 
micro-second  DBD  plasma  actuators.  Thermal  plasmas  are 
primarily  used  for  flow  control  from  transonic  to  super¬ 
sonic  to  hypersonic  flow  regimes,  and  include  counter-flow 
plasma  jets,  energy  deposition  methods  using  high  energy 
beams  like  electron  and  micro-wave,  and  gas  heating  using 
arc  or  electric  discharge  between  electrodes  at  the  vehicle 
surface.  Thermal  plasmas  have  been  reported  to  provide 
net  drag  and  surface  heat  transfer  reduction  by  affecting 
the  shock  thickness  and  speed.  Both  the  widening  of  the 
shock  and  the  increase  of  its  standoff  distance  from  the 
vehicle,  results  in  the  weakening  of  the  shock  intensity, 
thus  reducing  the  drag  and  surface  heat  transfer  to  the  vehi¬ 
cle.  It  is  a  requirement  for  such  mechanisms  to  have  energy 
efficiency  above  unity.  Despite  reported  success,  the 
energy  budget  of  this  approach  remains  high  (few  KWs  or 
higher ^^).  Though  the  cited  time  averaged  power  require¬ 
ment  maybe  few  hundred  watts,^"^  the  instantaneous  power 
requirement  may  exceed  few  kilowatts^’^’^^  to  even 

1 0  27 

megawatts.  ’ 

DBD  plasma  actuators  have  been  implemented  in 
nanosecond  and  microsecond  pulse  widths.  Nanosecond 
pulsed  DBD  plasma  actuators  have  also  been  successful  in 
high  speed  flow  manipulation.  As  an  example,  a  series  of 
high  voltage  nanosecond  pulses  producing  rapid  localized 
heating  are  used  to  reattach  the  flow  for  Mach  0.85  in  Ref. 
28.  In  quiescent  flow,  this  nanosecond  pulsed  DBD  plasma 
actuator  is  also  known  to  produce  shock  waves.  These 
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compression  waves  result  in  a  20%  increase  in  shock 
stand-off  distance  of  a  Mach  5  flow  over  a  cylinder/^ 
However,  the  advantage  of  heat  loss  due  to  the  increased 
shock  distance  is  seen  to  be  compensated  by  the  net  heat 
added.  In  fact,  the  profiles  for  Stanton  number,  Ch, 
increased  by  two  orders  of  magnitude.  Therefore,  nanosec¬ 
ond  pulsed  actuation,  though  promising  for  increasing  the 
shock  stand-off  distance,  seems  to  add  more  heat  during  its 
duty  cycle.  This  is  detrimental  to  the  re-entry  vehicle.  As 
an  alternative,  we  aim  to  explore  the  effect  of  micro¬ 
second  pulsed  DBD  plasma  to  re-entry  vehicle,  due  to  its 
low  power  requirement. 

According  to  Ref.  29,  micro-second  pulsed  DBD 
plasma  actuators  are  most  effective  at  bifurcation  points  in 
the  flow  field,  to  excite  the  instability  modes.  Their  applica¬ 
tions  are  limited  to  the  flow  velocities  up  to  60  m/s.  For 
hypersonic  flow  over  a  cylinder,  the  flow  velocities  are  very 
low  close  to  the  stagnation  point.  Thus,  it  allows  us  to  use 
micro-second  DBD  plasma  actuators,  where  the  main 
mechanism  is  momentum  addition  through  a  plasma  body 
force^^  to  a  flow  at  low  velocity  (<=50  m/s).  This  signifi¬ 
cantly  affects  the  net  heat  transfer  to  the  surface  of  the 
re-entry  vehicle  at  the  stagnation  point  by  locally  manipu¬ 
lating  the  flow.  In  addition,  the  power  requirement  for  this 
device  is  also  low,  making  the  powering  scheme  less 
complicated. 

This  paper  is  organized  as  follows.  Section  II  provides 
the  governing  equations  for  Navier  Stokes  equations  for  a 
single  species  and  multiple  species  relevant  for  hypersonic 
flows  with  and  without  full  thermo- chemical  non-equilib¬ 
rium.  Section  III  discusses  the  standard  plasma  DBD  actua¬ 
tor  and  the  formulation  of  the  body  force.  Section  IV 
describes  numerical  methodology  of  our  in-house  discontin¬ 
uous  Galerkin  (DG)  based  multi- scale  ionized  gas  (MIG) 
flow  code,  used  for  simulating  the  problems  presented  in 
this  paper.  Problem  description  and  boundary  conditions  for 
hypersonic  flow  problem  over  a  cylinder,  both  with  and 
without  thermo-chemistry,  is  presented  in  Sec.  V.  Section 
VI  presents  results  for  hypersonic  flows  without  thermo¬ 
chemistry  for  three  configurations  used  for  plasma  actua¬ 
tors.  Section  VII  summarizes  the  findings  for  hypersonic 
flows  with  full  non-equilibrium  flow  chemistry.  The  base¬ 
line  cases  are  compared  to  other  published  numerical 
results.  Based  on  these  studies,  the  most  receptive  actuator 
configuration  for  thermal  heat  flux  control  is  documented. 
Finally,  conclusions  are  made  about  the  key  findings  and 
future  directions. 


II.  GOVERNING  EQUATIONS 

We  consider  hypersonic  flows  both  with  and  without 
thermo-chemistry.  The  governing  equations  can  be  written 
as 


dU 

^  +  v.F.  =  v. 


(1) 


Here,  the  solution  variable,  U  and  the  inviscid  and  vis¬ 
cous  flux  terms,  i.e.,  and  respectively,  are  given 


’  p  ' 

pu 

pv 

pu 

pu^  +  F 

puv 

u  = 

pv 

•  = 

5  ^  1 

puv 

•  F^  = 

5  ^  1 

pv^  +  F 

.PE_ 

u[pE  +  F) 

_v(pF  +  F)_ 

0 

0 

^xx 

II 

'^yx 

'^xy 

Tyy 

UXxx  +  V^xy  +  qx  _ 

_UTyx  +  +  qy  _ 

(2) 


where  p  is  the  gas  density,  u  and  v  are  x-  and  y-velocities, 
E  is  the  total  energy,  P  is  the  pressure,  and  t  and  q  are  shear 
stress  and  the  conduction  heat  flux.  For  hypersonic  flows 
with  thermo-chemistry,  we  use  multi-species  Navier  Stokes 
equation.  A  general  form  of  this  equation  is  written  as 


dU 

dt 


+  V.Fi 


=  V.Fv  +  G. 


(3) 


The  solution  vector,  U,  the  inviscid  and  viscous  flux 
terms.  Ft  and  F^  and  the  source  term,  G  (in  Eq.  (3))  are  given 
in  Eq.  (4).  is  the  species  density,  with  s=  1,  2,...,  5  for  5 
species  model  of  O2,  NO,  N,  O,  and  N2.  p  is  the  total  density, 
which  is  the  sum  of  all  species  density,  pu  is  the  momentum 
(consisting  of  both  x  and  y  components),  pet  and  pey  are  the 
total  and  vibrational  energy,  respectively,  pet  is  the  sum  of 
the  vibrational  internal  energy,  translational-rotational  inter¬ 
nal  energy,  kinetic  energy  and  the  heat  of  formation  of  the 
present  species  (Eq.  (5)).  ht  is  the  total  enthalpy,  given  by  the 
sum  of  specific  internal  energy  plus  pressure  (Eq.  (6)) 
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where  hy  is  equal  to  Oy,  since  there  is  no  pressure  term  in  the 
vibrational  energy  equation,  as  the  kinetic  energy  is  consid¬ 
ered  only  in  the  total  energy  equation.  Vs,  appearing  in  the 
viscous  flux  term  for  the  continuity  equation  for  5  species 
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densities,  is  the  diffusion  velocity  of  the  species  and  is 
assumed  to  follow  Pick’s  law: 


p,Vs  =  -pDsVcs.  (7) 

Shear  stress,  t,  and  heat  fluxes,  q  and  for 
translational-rotational  modes  and  vibrational-electronic 
modes,  respectively,  are  given  by 

T  =  p{Vu  +  uV)  —  ^ pV  •  w/, 
q  =  —kVT, 

q^  =  -kVTy.  (8) 

Multi-component  diffusion  coefficient,  D^,  is  found 
using  Lewis  number,  Le  of  1.4.  Viscosity,  p,  is  found  using 
Wilke’s  semi-empirical  mixing  rule,^^  and  species  viscosities 
are  found  using  Blottner’s  model.^^  Conductivities,  k  and  k^ 
are  determined  from  an  Eucken  relation.^^  The  terms, 
^s^t,sPs^s  and  ^s^v,sPs^s  in  Efl-  (4)  in  the  viscous  flux 
term,  denote  the  transport  of  total  and  vibrational  energies, 
respectively,  due  to  mass  diffusion  of  the  different  species 
present.  Finally,  CDs  is  the  production  rate  of  the  different 
species,  due  to  the  chemical  reactions  identified  in  Eq.  (9). 
The  chemistry  source  terms  are  based  on  a  model  proposed 
by  Ref.  33.  The  first  three  reactions  are  the  dissociation  reac¬ 
tions  of  N2,  O2,  and  NO.  Here,  M  is  the  third  body  that  can 
be  any  of  the  5  species  present,  colliding  with  the  dissociat¬ 
ing  molecule 


V2+M^V  +  V  +  M 

V  +  O+M 
V2  +  O  ^  VO+V 

NO^O  ^  02^N.  (9) 

The  last  two  reactions  in  Eq.  (9)  are  the  exchange  reac¬ 
tions,  resulting  in  the  oxidation  of  N2  and  NO  species,  and 
production  of  the  atomic  nitrogen.  CDsCy^s  and  Qtv  in  the 
source  term,  G  in  Eq.  (4)  are,  respectively,  the  production/ 
destruction  of  the  vibrational  energy  (due  to  production/ 
destruction  of  species)  and  the  energy  exchange  between 
translational-rotational  and  vibrational  modes.  Landau- 
Teller  model  is  used  for  Qtv»  where  the  Landau-Teller 
inter-species  relaxation  time  is  given  in  Ref.  34.  The  energy 
first  transfers  from  kinetic  to  the  internal  translational- 
rotational  modes,  causing  a  rise  in  their  temperature  as  the 
flow  stagnates  in  a  hypersonic  flow.  The  internal  vibrational 
modes  will  then  equilibrate  with  the  other  modes  based  on 
the  modeled  time  scales  that  are  a  function  of  temperature 
and  pressure.  Further  relevant  details  about  this  model  can 
found  from  Ref.  35. 

III.  MICRO-SECOND  PULSED  DBD  PLASMA 
ACTUATORS 

DBD  plasma  actuators  have  been  used  for  many  flow 
control  applications  at  subsonic  Mach  number The 
actuator  consists  of  two  horizontally  displaced  electrodes 


separated  by  a  dielectric  material.  The  powered  electrode 
is  typically  exposed  to  the  surrounding  gas,  and  the 
grounded  electrode  embedded  beneath  the  dielectric. 
Application  of  high  sinusoidal  voltage  with  a  frequency 
of  1  50  kHz  and  amplitude  of  5  60  kVpp,  results  in  high 
electric  field  near  the  dielectric  surface.  This  high  electric 
field  ionizes  the  gas  next  to  the  electrodes,  and  the  combi¬ 
nation  of  electric  field  and  charged  particles  induces  body 
force  on  the  fluid  next  to  the  wall,  generating  a  wall  jet 
(Fig.  1). 

A  linear  configuration  is  the  standard  geometry  of  a 
DBD  plasma  actuator  (see  schematic  in  Fig.  1).  The  body 
force  distribution  for  the  linear  geometry  plasma  actuator 
was  modeled  in  Ref.  20,  using  first  principles  based  plasma 
discharge  simulations.  The  force  distribution  is  given  in 
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Here,  =  2.6,  Fy,o  =  2.0,  =  7.2  x  10\  Py  =  9x  10^ 

and  yb  =  0.00333.  Xq  and  yo  represent  the  edge  location  of  the 
electrodes.  Though  this  body  force  is  not  as  accurate  as  a  fully 
coupled  simulation  of  the  plasma  and  fluid  dynamics,  it  is 
more  accurate  than  phenomenological  models.  This  approach 
also  avoids  the  expensive  computations  required  for  fully 
coupled  simulations  that  include  the  electro-magnetic  field. 

The  force  distribution  for  the  plasma  actuator  imposed 
on  the  cylinder,  following  Eq.  (10),  is  shown  in  Fig.  2. 
Herein,  only  F^  is  shown  (here,  x  is  parallel  to  DBD  actuator, 
and  is  different  from  x  and  y  axis  shown  in  the  figure).  This 
Ex  is  multiplied  by  a  factor  of  15294  to  get  a  representative 
integrated  value  of  40  mN/m,  the  average  body  force 
obtained  from  plasma  DBD  actuators.  The  force  distribution 
is  on  the  flow  field  in  front  of  the  stagnation  region,  which  is 
exposed  to  an  incoming  Mach  number  of  17.65,  in  this  pa¬ 
per.  The  actuator  applied  to  this  case  is  also  shown  in  Fig.  2, 
and  it  actuates  the  fluid,  next  to  the  stagnation  point,  in  the 
upward  or  +y-direction. 


FIG.  1 .  Schematic  of  DBD  plasma  actuator,  with  exposed  electrode  (A)  and 
grounded  electrode  (C).  A  sinusoidal  AC  pulse  of  high  potential  difference 
across  the  exposed  and  grounded  electrodes  causes  a  high  electric  field  near 
the  surface,  which  causes  the  air  to  ionize.  This  electric  field  then  applies  a 
body  force  on  charged  particles  near  the  surface  that  causes  the  neutral  gas 
to  the  left  of  exposed  electrode  to  be  pulled  in  an  induced  wall  jet. 
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FIG.  2.  Force  distribution  of  Fx  (acting  in  the  positive  y  direction)  for 
plasma  actuator  applied  at  the  stagnation  point  of  the  cylinder. 


As  can  be  seen  in  this  figure,  the  plasma  body  force  is 
localized  at  high  magnitude.  Its  predominant  effect  is  con¬ 
centrated  within  a  region  of  3  10  mm  about  the  stagnation 
point.  Since  the  flow  velocities  are  very  low,  next  to  the  stag¬ 
nation  point,  a  significant  effect  of  the  plasma  actuator  on 
the  hypersonic  flow  is  to  be  expected  near  the  stagnation 
region. 

The  applicability  of  the  force  model  (which  is  generally 
applied  at  STP  conditions),  given  by  Eq.  (10),  at  hypersonic 
flow  conditions  at  the  stagnation  point  of  the  cylinder  is  an 
engineering  assumption,  and  can  be  justified  in  the  following 
manner.  The  flow  velocities  at  the  region  where  plasma  body 
force  acts  are  lower  than  15  m/s  (see  Fig.  14).  This  velocity 
regime  agrees  well  with  the  literature  of  DBD  plasma  flow 
control.  The  free  stream  pressure  for  the  given  flow  is  56  Pa, 
corresponding  to  a  50  km  re-entry  altitude.  Corresponding 
pressure  at  the  stagnation  point  is  0.227  atm.  Based  on  pub¬ 
lished  report,  a  stable  and  homogenous  discharge  can  be 
obtained  for  pressures  between  lO"^  and  10^  Pa,"^^  making  this 
application  feasible  at  the  given  pressure  at  the  stagnation 
point.  Additionally,  low  pressure  thrust  measurements"^^ 
also  show  similar  values  at  0.227  atm  and  1  atm.  This  justi¬ 
fies  the  use  of  the  force  model  for  the  given  pressure  condi¬ 
tions  at  the  stagnation  point.  Finally,  the  temperature  values 
rise  from  500  K  (at  the  wall)  to  nearly  10750  K  and  5250  K 
(at  about  5  mm  distance  away  from  the  cylinder  wall,  where 
the  plasma  body  force  is  significant,  see  Figs.  10  and  22, 
respectively),  for  hypersonic  flows  with  and  without  chemis¬ 
try,  respectively.  Even  though  the  temperature  value  is  far 
from  STP,  the  presence  of  high  energy  molecules  and  free 
radicals  (N2  and  O2  for  the  case  without  chemistry  and  N,  O, 
NO,  N2,  and  O2  for  the  case  with  chemistry),  near  the  cylin¬ 
der  surface,  will  only  aid  the  ionization  process.  Hence,  a 
better  performance  for  the  surface  heat  reduction  is  to  be 
expected  than  currently  predicted  using  the  plasma  DBD 
force  model  available  in  hand.  In  addition,  to  ascertain  its 
general  applicability,  this  force  model  has  been  validated 
with  physical  experiments  under  quiescent  and  subsonic  flow 
conditions. 


IV.  NUMERICAL  METHODOLOGY 

Our  in-house  MIG  solver  was  used  for  simulating  the 
problems  presented  in  this  work.  MIG  is  a  multi-physics 
finite  element  based  code.  It  has  been  used  to  solve  various 
fluid-plasma  interaction  problems. MiG’s  original 
formulation  is  Galerkin  based  finite  element  method,  which 
requires  the  solution  to  be  continuous  across  the  elements’ 
interfaces  (Fig.  3).  This  however  limits  the  application  of 
MIG  to  diffusion  dominated  problems.  For  extending  MIG 
to  high-speed  flows,  a  DG  framework  was  incorporated.  This 
was  conveniently  done  owing  to  MiG’s  modular  structure. 
Some  successful  applications  of  DG  method  applied  to  MIG 
have  been  presented  in  Refs.  47,  48,  and  35. 

The  discontinuous  Galerkin  method  is  a  finite  element 
based  method,  where  in  the  solution  continuity  constraint 
across  elements’  interfaces  is  relaxed.  This  makes  DG 
method  suitable  for  various  high-speed  flow  problems. 

DG  methods  are  known  to  combine  the  relative  advantages 
of  finite  volume  method  (for  convection  dominated  prob¬ 
lems)  and  finite  element  methods  (for  unstructured  meshes 
and  complex  geometries).  In  addition,  DG  method  is  also  a 
compact  stencil  high-order  accurate  method  that  makes  it 
convenient  to  extend  to  higher  orders  without  any  extra 
coding  effort. 

For  basis  functions,  we  use  2D  Jacobi  polynomials^^ 
given  by 


x{(l-,)-+"(l +,)"+■}.  (11) 

To  summarize  the  DG  method  used  in  this  work,  local 
Fax-Friedrichs  flux^^  and  second  formulation  of  Bassi  and 
Rebay,  namely,  BR2  scheme^^  are  used  for  the  inviscid  and 
viscous  numerical  fluxes,  which  are  defined  at  the  element 
interfaces.  Backward  Euler,  an  implicit  time  integration 
method,  is  used  for  advancement  of  solution  with  high 
Courant-Friedrichs-Lewy,  CFL  numbers  (approximately 
5000). 

To  capture  the  shock,  we  use  r-p  adaptivity  method  pro¬ 
posed  and  presented  in  Refs.  35  and  48.  In  this  approach,  we 
use  pressure  based  shock  sensors^"^  (Eq.  (12))  to  clearly  iden¬ 
tify  the  shock  location 


CG 


DG 


FIG.  3.  Schematic  showing  solution  variation  across  two  neighboring  ele 
ments  in  continuous  Galerkin  (CG)  and  DG.  Solution  across  neighboring 
elements  is  required  to  be  continuous  in  CG  and  allowed  to  be  discontinuous 
in  DG. 
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FIG.  4.  Mesh  resolution  using  r  p 
adaptivity  iterations,  (a)  Mesh  before 
r  p  adaptivity.  Elements  in  the  shock 
are  identified  in  red  and  element  out 
side  the  shock  are  identified  in  blue, 
(b)  Mesh  after  first  r  p  adaptivity  itera 
tion.  Elements  outside  the  shock  are 
moved  to  within  the  shock  to  refine  the 
elements  in  the  shock,  (c)  Second  r  p 
adaptivity  iteration  and  (d)  third  r  p 
adaptivity  iteration.  Shock  mesh  reso 
lution  is  successively  improved  with 
each  iteration  of  r  p  adaptivity,  without 
increasing  the  total  number  of 
elements. 


Sk  —  log  10 


also  implemented  parallelization  of  MIG  using  two  parallel 
(12)  solver  platforms,  HYPRE^^  and  PETSc.^^  These  simulations 
are  run  using  PETSc  parallel  solver. 


where  [Ph]  in  Eq.  (12)  denotes  the  jump  in  the  pressure 
across  the  edges  of  an  element,  k,  and  {Ph}  denotes  the  aver¬ 
age  of  pressure  across  the  edge,  9Qk.  Sk  is  evaluated  for  all 
the  edges  of  an  element,  Qk  in  the  domain.  Eollowing  criteria 
(Eq.  (13))  is  used  to  identify  the  regions  of  the  shock  from 
the  smooth  flow  regions.  Shock  locations  are  identified  for  Sk 
greater  than  a  threshold  value  of  0.3  (this  gives  optimum  per¬ 
formance)  and  smooth  flow  regions  are  identified  for  <  0.3 

>  0.3;p  =  0,  shock_region 
Sk  <  0.3; p  >  1 ,  smooth_region. 

A  switch  “lim”  is  used  to  identify  shock  regions  and  an 
example  is  shown  in  Eig.  4  for  hypersonic  flow  around  a  cyl¬ 
inder.  p  =  0  is  the  red  region  indicating  shock  and  p  >  1  is 
used  in  the  blue  region,  indicating  smooth  flow  regions. 
Solution  accuracy  in  smooth  flow  regions  can  be  improved 
by  increasing  the  order  of  polynomial  approximation,  p, 
without  changing  the  mesh  resolution  there.  This  is  the 
advantage  afforded  by  DG  method.  But,  the  shock  region  is 
limited  to  p  =  0  and  hence  first  order  accurate.  A  reasonable 
way  to  improve  solution  accuracy  in  the  shock  region  is  to 
refine  the  mesh.  One  approach  taken  in  Refs.  54  and  55  is 
applying  h-refinement  in  the  shock  regions.  This  improves 
the  solution  accuracy  in  the  shock  regions.  However,  this 
also  results  in  large  number  of  elements  in  the  region  around 
the  shock,  which  may  not  lie  in  the  shock.  Eor  this  purpose, 
we  use  the  r-adaptivity  approach,  where  we  cluster  the  ele¬ 
ments  in  the  domain  outside  of  the  shock,  to  a  region  very 
close  to  the  shock.  This  refinement  is  based  on  the  shock 
location,  and  hence  applicable  to  all  body  shapes.  Successive 
r-p  adaptivity  iterations  resulting  in  increasingly  finer  mesh 
resolution  at  the  shock  are  shown  in  Eig.  4(b)  4(d).  We  have 


V.  PROBLEM  DESCRIPTION 

In  this  work,  we  consider  the  effect  of  plasma  DBD 
actuators  on  a  Mach  17  hypersonic  flow  over  a  cylinder,  of 
radius  1  m,  both  with  and  without  thermo-chemistry  (pre¬ 
sented  in  Secs.  VI  and  VII,  respectively).  Air  is  considered 
as  a  perfect  gas  for  viscous  hypersonic  flow,  and  a  mixture 
of  N2  and  O2  for  non-equilibrium  hypersonic  flow.  Eor  both 
simulations,  with  and  without  thermo-chemistry,  a  free- 
stream  temperature  of  200  K  and  a  fixed  cylinder  wall  tem¬ 
perature  of  500  K  is  assumed.  The  free-stream  and  boundary 
conditions  are  summarized  in  Table  I,  for  flows  with  and 
without  thermo-chemistry. 

Eor  the  non-equilibrium  flow,  the  mass  fractions  of  N2 
and  O2  are  in  a  ratio  of  0.7562:0.2438  in  the  free-stream.  Air 
is  also  considered  to  be  in  thermo-chemical  equilibrium  in 
the  free-stream,  hence  the  vibrational  temperature  for  air  is 
equal  to  its  translational-rotational  temperature,  i.e., 
Tv  =  200K.  Super-catalytic  wall  forces  the  flow  species, 
near  the  wall,  into  thermo-chemical  equilibrium.  Hence,  the 
vibrational  temperature  at  the  wall  is  equal  to  the  fixed  wall 


TABLE  I.  Free  stream  and  boundary  conditions  for  Ml 7  hypersonic  flow 
over  cylinder. 


Free  stream  and 

boundary  conditions  Without  thermo  chemistry  With  thermo  chemistry 

Mach  number,  M^o 

17.65 

17.65 

Density,  (kgM3) 

0.001 

0.001 

X  velocity,  u^o  (m/s) 

5000 

5000 

Toe  (K) 

200 

200 

Twall  (K) 

500  (fixed 

500  (fixed  temperature 

temperature  wall) 

and  super  catalytic  wall) 

u,  V  ((2)  wall)  (m/s) 

0  (no  slip  wall) 

0  (no  slip  wall) 
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1.? 
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FIG.  5.  Boundary  conditions  summarized  for  the  four  edges  of  the  domain. 
Free  stream  inlet  is  at  the  left  of  the  domain  and  the  cylinder  wall  at  the 
right.  Outflow  boundary  conditions  are  used  on  the  two  edges  on  top  and 
bottom  of  cylinder. 

temperature  of  500  K,  and  all  the  species  recombine  at  the 
wall  to  give  free-stream  species  at  the  wall. 

The  boundary  conditions  are  shown  in  Fig.  5.  The  out 
flow  is  set  to  be  supersonic,  but  the  boundary  layer  is  not. 
Having  incorrect  boundary  layer  profile  for  the  outflow, 
however,  will  not  affect  results  near  the  stagnation  region. 

VI.  EFFECT  OF  PLASMA  DBD  ACTUATORS  ON 
VISCOUS  HYPERSONIC  FLOWS 

First,  we  look  into  the  effect  of  micro- second  pulsed 
DBD  plasma  actuators  on  a  viscous  perfect  gas  hypersonic 
flow.  The  body  force,  as  a  result  of  plasma  actuation,  is  first 
applied  at  the  stagnation  point  on  the  cylinder  surface.  We 
first  note  the  effect  of  this  actuation,  on  the  heat  transfer 
coefficient,  Ch,  in  order  to  enhance  our  understanding  of  the 
physics  of  this  problem.  Various  designs  for  plasma  actua¬ 
tors  are  then  considered,  by  varying  locations  and  polarity 
(i.e.,  direction  of  the  force)  of  the  plasma  actuation.  This 
grid  has  256  x  100  elements  in  the  domain,  where  the  256 
elements  are  uniformly  distributed  in  the  9  (or  circumferen¬ 
tial)  direction  and  100  elements,  with  varying  sizes,  are 


FIG.  6.  Mesh  used  for  plasma  actuation  of  viscous  hypersonic  flow  has 
256  X  100  elements.  This  is  achieved  after  4th  r  p  adaptivity  iteration.  This 
mesh  shows  very  fine  mesh  obtained  at  the  shock. 

distributed  along  the  radial  direction.  A  part  of  final  r-p 
adapted  grid  is  shown  in  Fig.  6. 

We  consider  three  configurations,  as  shown  in  Fig.  7,  for 
the  placement  of  plasma  DBD  actuators  on  the  cylinder. 
“Plasmal”  refers  to  the  first  configuration  with  the  plasma  ac¬ 
tuator  placed  at  the  stagnation  point,  with  body  force  in  the 
+y-direction.  For  the  second  and  third  configurations,  two 
plasma  actuators  are  placed  symmetrically  on  the  cylinder 
surface,  at  a  distance  of  25  mm  (measured  along  the  circum¬ 
ference  of  the  cylinder)  away  from  the  stagnation  point  in  the 
positive  and  negative  y-directions.  In  the  second  configura¬ 
tion  or  “plasma2,”  both  plasma  actuators  apply  body  force 
away  from  each  other,  i.e.,  plasma  actuator  in  the  positive 
part  of  y-axis,  applies  force  in  the  +y-direction  and  the 
plasma  actuator,  in  the  negative  part  of  y-axis,  applies  force 
in  -y-direction.  In  the  third  configuration,  i.e.,  “plasmaS,”  the 
directions  of  both  forces  are  reversed  (from  “plasma2”),  so 
that  both  the  forces  act  towards  each  other.  The  schematic  in 
Fig.  7  illustrates  the  three  configurations.  The  effect  of  these 
three  configurations  on  the  Ch  profile  is  summarized  in  Fig.  8. 

Ch  measures  the  ratio  of  heat  transferred  to  the  solid  sur¬ 
face  (from  the  fluid)  to  the  thermal  capacity  of  the  fluid.  Its 
formula  is  given  in 


Fb 


plasmal  plasma^  plasmaS 

FIG.  7.  Three  configurations  for  plasma  DBD  actuators.  In  the  first  configuration,  plasma  actuator  is  applied  at  the  stagnation  point,  with  flow  forcing  in  the 
positive  y  direction.  Second  and  third  configurations  have  same  placements  of  the  actuators,  i.e.,  ±25  mm  distance  from  the  stagnation  point  along  the  circum 
ference  of  the  cylinder,  but  with  reversed  flow  forcing.  Second  configuration  has  forces  pointing  away  from  each  other  and  third  configuration  has  forces  point 
ing  towards  each  other. 
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FIG.  8.  Effect  of  plasma  actuation  on  the  surface  heat  transfer  coefficient,  Ch  of  hypersonic  flow  over  cylinder  problem  shown  in  (a)  zoomed  in  section  and  in 
(b)  full  plot.  The  case  of  “no  plasma”  from  MIG,  compares  well  to  profile  predicted  from  LAURA.^^  The  comparisons  of  Ch  profile  for  “plasmal,” 
“plasma2,”  and  “plasmaS”  are  also  shown.  Best  reduction  of  Ch  profile  is  obtained  for  plasmaS. 


Ch  =  y^^,  (14) 

2  P 00^00 

where  (q)^  is  the  net  heat  flux  in  the  wall  normal  direction. 

The  plasmal  conflguration  leads  to  both  an  increase  and 
reduction  of  Ch  profile  along  the  surface  of  the  cylinder.  A 
maximum  value  of  0.00846,  occurs  at  0=— 0.0136  (very 
close  to  stagnation  point),  and  minimum  value  of  0.00758 
occurs  at  0  =  0.025.  Since  the  radius  is  1  m,  these  locations 
correspond  to  —13.6mm  and  25  mm  distance  away  from 
stagnation  point  (along  the  cylinder  circumference).  Upon 
performing  integration  over  the  entire  surface,  we  can  see 
that  the  plasma  actuator  in  the  “plasmal”  conflguration  pro¬ 
vides  net  cooling  effect.  The  other  two  configurations, 
plasma2  and  plasma3  are  similar  to  moving  the  plasma  DBD 
effect  on  Ch  profile  (in  plasmal)  in  ±25  mm  direction  to  the 
location  of  plasma  actuator,  and  reversing  the  Ch  profile 
based  on  the  direction  of  the  plasma  body  force.  Thus,  the 


second  and  third  configurations  give  overall  heat  transfer  in¬ 
crement  and  reduction,  respectively. 

To  understand  the  flow  physics,  the  effect  of  the  plasma 
actuator  on  the  temperature  and  velocity  profiles,  near  the  wall 
(in  direction  normal  to  the  wall),  is  analyzed,  w.r.t.  “no  plasma” 
case.  For  “plasmal”  conflguration,  the  effect  of  the  actuator  is 
shown  in  Fig.  9,  both  with  and  without  plasma.  Visually,  the 
temperature  contours  are  pushed  closer  to  the  cylinder  surface 
for  negative  y-axis  and  away  from  the  surface  for  positive 
y-axis,  for  “plasmal”  relative  to  “no  plasma”  case.  Thus,  the 
actuator  reduces  the  Ch  value  at  about  25  mm  distance  down¬ 
stream  from  its  actual  location,  by  strongly  pushing  the  temper¬ 
ature  contours  away  from  the  cylinder  surface. 

It  can  therefore  be  seen  that  the  micro-second  pulsed 
plasma  DBD  actuator,  which  are  generally  considered  only 
appropriate  for  subsonic  flows,  have  significant  effect  on 
hypersonic  flows  (near  the  boundary  layer),  in  terms  of  heat¬ 
ing.  This  is  due  to  the  large  stagnation  region  in  the  hyper¬ 
sonic  flow,  with  low  velocities  next  to  the  wall  that  allows 


x(m) 


T(K) 

12620 

12580 

12200 

11000 

9000 

7000 

5000 

3000 

1000 


1  t  i  ^  1  t  1 1  i  1 

-0.95  -0.9 


1.2  -1.15  -1.1  -1.05  -1 

x(m) 


T(K) 

12620 

12580 

12200 

11000 

9000 

7000 

5000 

3000 

1000 


-0.95  -0.9 


FIG.  9.  Effect  of  “no  plasma”  vs  plasma  actuator  in  “plasmal”  configuration  on  the  temperature  contours  close  to  the  stagnation  region  of  the  cylinder. 
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FIG.  10.  Temperature  plots  in  the  flow,  shown  with  distance  (in  the  flow)  away  from  the  wall,  denoted  as  x^.  (a)  Overall  plot  for  different  y  locations  (shown 
in  the  legend)  and  (b)  zoomed  in  section.  Y  locations  correspond  to  the  distances  traveled  along  the  cylinder  surface  in  the  positive  or  negative  direction  from 
the  stagnation  point.  Solid  lines  refer  to  the  results  without  plasma  and  dotted  lines  refer  to  the  results  with  plasma  (plasmal). 


the  plasma  actuator  to  make  some  significant  effect  on  the 
surface  heating  in  hypersonic  flow,  by  adding  momentum 
into  the  flow. 

Next,  we  look  at  the  temperature  and  velocity  profiles 
(Figs.  10  and  11),  near  the  cylinder  surface,  as  we  move 
away  from  the  wall  in  the  normal  direction.  We  notice  from 
Fig.  10  that  all  the  temperature  profiles  fall  into  a  single 
curve  for  the  “no  plasma”  case.  For  “plasmal”  case,  the  tem¬ 
perature  profiles  are  above  and  below  relative  to  the  profiles 
of  the  “no  plasma”  case  for  locations  closer  to  the  actuation 
point  and  in  the  positive  and  negative  half  of  y-axis,  respec¬ 
tively.  This  corresponds  to  the  increment  and  decrement 
observed  in  Ch  profile  in  Fig.  8.  From  Fig.  10,  we  note  that 
the  maximum  value  of  Ch  is  at  —12.5  mm,  close  to 


—  13.6  mm,  seen  in  Fig.  8.  Velocity  profiles  in  Fig.  11  indi¬ 
cate  significant  increase  in  y-velocity,  from  1  m/s  to  8.5  m/s, 
at  a  distance  of  nearly  1  mm  away  from  cylinder  surface  at 
the  stagnation  point. 

The  temperature  contours  for  second  and  third  configu¬ 
ration  (close  to  the  stagnation  region)  are  shown  in  Fig.  12. 
In  “plasma2”  configuration,  the  temperature  contours  are 
compressed  closer  to  the  surface,  relative  to  “no  plasma” 
case,  from  y  =  — 50  mm  to  +50  mm.  Thus,  Ch  value  for 
plasma2  configuration  is  higher  than  “no  plasma”  case  for 
0  =  — 0.05  0.05.  For  the  rest  of  the  region,  Ch  profile  falls 
below  the  Ch  profile  for  the  “no  plasma”  case. 

Correspondingly  for  “plasmaS”  configuration,  the  tem¬ 
perature  contours  are  pushed  away  from  the  cylinder  wall. 


FIG.  1 1.  X  and  Y  velocity  profiles  shown  with  x^,  the  distance  measured  in  the  wall  normal  direction,  away  from  the  wall.  Y  velocity  is  almost  parallel  to  the 
cylinder  surface,  near  the  stagnation  point,  and  is  thus  relevant  to  understand  the  cause  of  observed  change  in  Ch  value  near  stagnation  point. 
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FIG.  12.  Temperature  contours  near  stagnation  region  for  two  plasma  actuators  at  25  mm  (w.r.t.  stagnation  point  on  cylinder  surface),  acting  (a)  with  body 
forces  pointed  away  from  each  other,  i.e.,  plasma2  configuration  and  (b)  with  body  forces  pointed  towards  each  other,  i.e.,  plasmaS  configuration. 


for  y  =  —50  mm  to  +50  mm,  thus  causing  overall  reduction 
in  Ch  value  relative  to  “no  plasma”  case.  What  causes  the 
temperature  contours  to  be  pulled  into  or  pushed  away  from 
the  wall  in  plasma2  and  plasmaS  conhgurations  is  explained 
through  the  velocity  schematic  given  in  Fig.  13,  for 
“plasma  1”  conhguration.  Plasma  body  force  induces  a  flow 
(shown  in  red)  in  the  direction  of  the  force.  The  oncoming 
flow  (shown  in  blue)  hits  the  cylinder  surface  (approximated 
by  a  flat  wall  in  black).  Without  the  plasma  body  force,  the 
flow  has  zero  tangential  velocity  (y-velocity)  a  few  milli¬ 
meter  distance  away  from  the  wall,  and  hence  maximum 
value  of  Ch  is  experienced  at  the  stagnation  region.  The  flow 
induced,  by  the  plasma,  causes  a  leftward  shift  (from  the 
stagnation  point)  of  the  point  of  zero  tangential  velocity. 
This  causes  the  maximum  in  Ch  to  shift  slightly  to  the  left  in 
Fig.  8.  The  question,  that  arises,  is  why  Ch  profile  shoots 
above  the  maximum  value  of  Ch  at  the  stagnation  point  for 
“no  plasma”  case. 

Velocity  profiles,  induced  by  the  plasma  actuator  in  the 
“plasma  1”  configuration,  relative  to  “no  plasma”  case  are 
shown  in  Fig.  14.  Actuator  decelerates  the  flow  (going  in  the 
negative  direction),  to  the  left  of  the  stagnation  point,  hence 
causing  increased  heating  at  the  surface  there,  and  acceler¬ 
ates  the  flow,  to  the  right  of  the  stagnation  point,  hence 


Flow  directiori 

Flow  induct  by  plasma 


FIG.  13.  Schematic  of  flow  induced  (in  red)  by  plasma  body  force  with  force 
direction  in  direction.  The  oncoming  flow  (in  blue)  hits  the  wall  (in 
black).  The  cylinder  surface  close  to  the  stagnation  region  can  be  approxi 
mated  by  a  flat  wall. 


reducing  heat  transfer  to  that  surface.  Suction  of  the  high 
temperature  gas  causes  increased  heat  transfer,  below  the 
stagnation  point.  Added  momentum  pushes  away  the  high 
temperature  gas  (from  the  wall),  resulting  in  reduced  heat 
transfer,  above  the  stagnation  point.  Temperature  profiles  for 
plasma2  and  plasmaS  configurations  in  comparison  to  “no 
plasma”  case  are  given  in  Fig.  15.  It  can  be  seen  above  that 
temperature  profiles  for  the  “plasma2”  case,  primarily  lie 
below  the  temperature  profile  for  “no  plasma”  case  and  that 
for  “plasmaS”  case,  primarily  lie  above  the  “no  plasma” 
case,  thus  being  consistent  with  higher  and  lower  Ch  profiles 
observed  for  “plasma2”  and  “plasmaS”  cases,  respectively, 
in  reference  to  “no  plasma”  case. 

Same  observation  is  also  made  in  the  y-velocity  profiles 
for  “plasma2”  and  “plasmaS”  cases  (see  Fig.  16).  Here,  the 
plasma  DBD  actuator  located  in  the  positive  y-axis. 


FIG.  14.  Y  velocity,  v  (m/s)  (parallel  to  the  cylinder  wall)  is  generated  for 
several  locations  close  to  the  stagnation  point,  lying  both  to  the  left  and  the 
right.  This  is  plotted  for  both  with  and  without  plasma  actuators  (solid  lines 
represents  the  case  without  plasma  actuator  and  dotted  lines,  the  case  with 
the  plasma  actuator).  This  plot  is  for  comparison  of  “plasmal”  case  to  “no 
plasma”  case. 
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FIG.  15.  Temperature  profiles  for  “plasma2”  and  “plasmaS”  configurations  in  comparison  to  the  “no  plasma”  case.  For  plasma2  case,  the  temperature  profiles 
(shown  in  dotted  profile)  have  steeper  gradient  than  “no  plasma”  case  and  for  “plasmaS”  case,  the  temperature  gradient  has  lower  gradient. 


introduces  a  wall  jet  in  the  direction  of  the  fluid  flow  in  the 
“plasma2”  configuration  and  in  the  opposite  direction  of  the 
fluid  flow  in  the  “plasmaS”  configuration.  By  increasing  the 
flow  velocity  next  to  the  wall,  the  plasma  actuator  is  able  to 
entrain  more  of  the  high  temperature  fluid,  thus  producing 
more  surface  heating  for  “plasma2”  configuration,  and  by 
decreasing  the  flow  velocity  next  to  the  wall,  it  is  able  to 
push  the  high  temperature  fluid  away  from  the  wall,  thus 
reducing  heating  next  to  the  wall  for  plasmaS  configuration. 

Vll.  EFFECT  OF  PLASMA  DBD  ACTUATORS  ON 
NON-EQUILIBRIUM  HYPERSONIC  FLOWS 

After  successfully  looking  into  the  effect  of  plasma  DBD 
actuator  on  the  Mach  17  flow  over  hypersonic  flow,  we  inves¬ 
tigate  its  effect  on  a  cylinder  in  hypersonic  flow  with  thermo¬ 


chemical  non-equilibrium.  But  before  this,  we  show  basic 
comparison  of  results  for  MIG  with  the  US3D  code  devel¬ 
oped  at  the  University  of  Minnesota.^^  US3D  is  an  unstruc¬ 
tured  CFD  code  for  hypersonic  flow  solution  used 
extensively  to  solve  many  different  hypersonic  flow  problems 
and  has  gone  through  extensive  validation  and  verification. 
The  mesh  used  in  MIG,  for  thermo-chemical  non-equilibrium 
hypersonic  flow  simulation,  is  shown  in  Fig.  17.  This  was 
obtained  after  the  4th  r-p  adaptivity  iteration. 

Temperature  contours  for  N2  +  O2  case  are  shown  in 
Fig.  17  for  both  MIG  and  US3D.  Both  the  results  compare 
very  well  to  each  other.  The  maximum  temperature  in  MIG 
is  found  to  be  7700  K  and  in  US3D  to  be  6500  K.  This  differ¬ 
ence  results  from  finer  mesh  used  in  MIG,  from  the  succes¬ 
sive  (4th)  r-p  adaptivity  iteration,  in  comparison  to  coarser 
mesh  being  used  in  US3D.  The  temperature  contours  show  a 


FIG.  16.  Y  velocity  profiles  for  “plasma2”  and  “plasmaS”  configurations  in  comparison  to  “no  plasma”  case.  In  plasma2  configuration,  flow  is  accelerated, 
and  in  plasmaS  configuration,  flow  is  decelerated  relative  to  “no  plasma”  case. 
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FIG.  17.  Mesh  and  comparisons  of 
temperature  contours  for  Mach  17  non 
equilibrium  hypersonic  flow,  (a)  Mesh 
obtained  after  4th  r  p  adaptivity  itera 
tion  for  MIG,  (b)  temperature  contours 
comparison  obtained  from  MIG  code, 
and  (c)  US3D  code.  Both  results  com 
pare  very  well  to  each  other. 
Maximum  temperature  in  MIG  is 
7700  K  and  in  US3D  is  6500  K,  and 
the  difference  is  due  to  finer  mesh  at 
the  shock  in  MIG  obtained  with  sue 
cessive  (4th)  r  p  adaptivity  iteration  in 
comparison  to  coarser  mesh  at  the 
shock  for  US3D. 


spiked  nature  in  the  shock  and  a  smoothed  out  nature  to  the 
inside  of  the  shock  layer.  This  results  from  strong  thermo¬ 
chemical  non-equilibrium  within  the  shock,  which  quickly 
equilibrates  in  the  shock  layer,  due  to  large  number  of  mo¬ 
lecular  collisions. 

Contours  for  vibrational  temperature  also  compare  very 
well  for  MIG  and  US3D  (see  Fig.  18).  For  Ty,  the  maximum 
temperature  in  flow  field  for  MIG  is  6300  K  and  for  US3D  is 
6550  K.  Hence,  the  maximum  values  of  vibrational  tempera¬ 
ture  are  much  closer  for  both  the  codes  than  it  is  for  the  maxi¬ 
mum  values  of  translational-rotational  temperature.  Contours 


for  mass  fraction  of  O2  are  shown  in  Fig.  18,  and  both  codes 
show  complete  dissociation  of  O2  in  the  shock  layer.  This  is 
due  to  temperatures  exceeding  2000  K.  Corresponding  results 
for  mass  fractions  of  O  are  seen  in  Fig.  19.  Dissociation  of  N2 
into  N  atoms  is  incomplete  with  maximum  dissociation  (of 
nearly  0.07  of  total  mass  fraction,  see  Fig.  19)  occurring  just 
after  the  shock.  Significant  dissociation  of  N2  requires  higher 
temperatures,  of  the  order  of  7000  K  or  higher.  A  small  por¬ 
tion  of  the  atomic  nitrogen,  N  (produced  by  the  dissociation), 
combines  with  atomic  oxygen,  O  to  form  NO  close  to  the 
shock  region  (Fig.  20).  Away  from  the  shock,  in  the  shock 


FIG.  18.  Contours  for  vibrational  temperature  for  (a)  MIG  and  (b)  US3D,  and  mass  fraction  of  O2  for  (c)  MIG  and  (d)  US3D.  Both  results  compare  very  well 
to  each  other.  Maximum  vibrational  temperature  in  MIG  is  6300  K  and  in  US3D  is  6550  K,  and  the  difference  is  again  attributed  to  meshing  difference  in  both 
the  codes.  For  O2,  a  complete  dissociation  is  noticed  in  both  the  codes. 
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FIG.  19.  Mass  contours  of  O  for  (a)  MIG,  (b)  US3D  and  N2  for  (c)  MIG  and  (d)  US3D.  Comparisons  are  fairly  accurate.  Maximum  mass  fraction  for  O,  i.e., 
0.2438,  is  found  close  to  the  top  and  bottom  regions  of  the  cylinder.  This  correlates  to  the  mass  fraction  of  NO,  in  that  the  maximum  O  is  observed,  when  NO 
(formed  in  the  shock)  dissociates  into  N  and  O  (due  to  lower  temperature  in  the  shock  layer)  giving  extra  O.  A  partial  dissociation  of  N2  is  observed  close  to 
the  stagnation  region,  as  the  temperature  is  not  of  order  of  7000  K  or  higher  for  significant  dissociation. 


layer,  the  temperatures  are  not  high  enough  to  sustain  NO 
which  again  dissociates  into  atomic  N  and  O.  Thus,  in  the 
contours  of  N,  we  see  the  maximum  value  of  mass  fraction 
being  just  after  the  shock  rather  than  in  the  shock. 

As  N2  and  O2  continue  to  dissociate  past  the  shock,  the 
temperature  contours  which  have  a  maximum  in  the  shock, 
continue  to  drop  unlike  purely  viscous  case,  where  the 


temperature  increases  as  the  flow  stagnates  until  the  viscous 
boundary  layer  is  reached.  As  we  go  away  from  the  stagna¬ 
tion  line  towards  either  the  top  or  bottom  of  cylinder,  the 
temperature  values  decrease,  and  hence  the  amount  of  NO 
(dissociating  into  N  and  O)  also  decreases,  thus  we  see  that 
maximum  mass  fraction  of  O  is  at  top  and  bottom  regions  of 
cylinder  and  not  near  the  stagnation  line. 
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FIG.  20.  Mass  fraction  contours  of  NO  for  (a)  MIG  and  (b)  US3D  and  of  N  for  (c)  MIG  and  (d)  US3D.  Maximum  mass  fraction  for  NO  and  its  production 
occurs  in  the  shock  by  combination  of  atomic  N  and  O  produced  by  dissociation  reactions  of  N2  and  O2  molecules.  As  the  temperatures,  being  maximum  at  the 
shock,  drop  down  in  the  shock  layer,  NO  again  dissociates  back  into  N  and  O  atoms.  Production  of  atomic  N  corresponds  to  the  dissociation  of  N2  shown  in 
Fig.  19. 
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FIG.  21.  The  effect  of  the  plasma  DBD  actuator  on  Ch  profile,  for  non 
equilibrium  hypersonic  flow  is  shown.  In  reference  to  the  base  case,  with  no 
plasma  (shown  in  red),  plasma  1  case  gives  both  increase  and  decrease  in  Ch 
value,  plasma2  case  gives  overall  increase  in  Ch  value  and  plasmaS,  as 
expected,  gives  overall  reduction  of  heating. 

Having  verified  our  non-equilibrium  hypersonic  fiow 
code  with  US3D,  we  now  investigate  the  effect  of  the  plasma 
actuators  (in  three  configurations  shown  in  Fig.  7)  on  the 
non-equilibrium  hypersonic  fiow.  Just  as  in  Sec.  VI,  we  look 
at  the  effect  on  the  surface  heating  profile  (Fig.  21).  The 
effect  is  similar  to  Fig.  8,  which  is  expected.  In  non¬ 
equilibrium  hypersonic  fiow,  the  heat  fiux  into  the  wall  is 


equal  to  net  conduction  based  on  both  translational- 
rotational  and  vibrational  temperatures,  i.e.,  T  and  T^,  and 
the  enthalpy  transport  of  diffusing  species  close  to  the  wall. 
Comparing  plasmaS  configuration  to  the  base  case  of  no 
plasma,  reduction  of  temperature  gradient  (at  stagnation 
line)  next  to  the  wall  is  observed  (Fig.  22).  This  is  true  for 
both  translational-rotational  temperature,  T,  and  vibrational 
temperature,  T^.  The  effect  is  dominant  to  within  10  mm  dis¬ 
tance  away  from  the  wall.  In  addition  to  the  temperature  pro¬ 
files,  the  application  of  plasma  DBD  actuator  is  also  seen  to 
have  effect  on  the  species  mass  fractions.  This  effect  is  not 
so  pronounced  for  O2  and  O,  as  can  be  seen  in  Fig.  23. 

However,  a  more  significant  effect  is  seen  on  stagnation 
line  plots  for  N  and  NO  (in  Fig.  24).  The  difference  of 
plasma3  case,  w.r.t.  the  base  case  is  seen  up  to  25  mm  dis¬ 
tance  away  from  the  wall.  In  Fig.  24(a),  we  notice  sudden 
reduction  of  N,  close  to  the  wall,  for  both  with  and  without 
the  plasma  actuator.  This  happens  from  a  reduction  in  tem¬ 
perature,  close  to  the  wall  (due  to  the  boundary  condition). 
At  lower  temperatures,  nitrogen  cannot  sustain  in  the  atomic 
form,  and  hence  recombines  to  N2  (primarily)  and  NO  near 
the  wall.  This  recombination  releases  large  amount  of  heat¬ 
ing  to  the  wall.  A  minor  change  in  temperature,  because  of 
plasma,  causes  significant  variation  in  mass  fraction  of  N,  up 
to  distance  of  25  mm  away  from  wall.  For  NO,  there  is 
increase  in  mass  fraction,  close  to  the  wall  (from  4  mm  to 
2  mm  next  to  the  wall).  This  is  due  to  atomic  N  and  O 


FIG.  22.  Effect  of  plasmaS  configura 
tion  on  stagnation  line  plots  of  T  and 
Tv  is  shown  next  to  the  wall,  in  refer 
ence  to  “no  plasma”  case.  Only  a 
minor  effect  is  seen  in  the  temperature 
plots  for  distance  up  to  10  mm  away 
from  cylinder  wall. 


FIG.  23.  Effect  of  plasmaS  configura 
tion  on  stagnation  line  plots  of  O2  and 
O,  is  shown  next  to  the  wall,  in  refer 
ence  to  “no  plasma”  case.  Only  a 
minor  effect  is  seen  in  the  plots  for  dis 
tance  up  to  6  mm  away  from  cylinder 
wall. 
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FIG.  24.  Effect  of  plasmaS  configura 
tion  on  stagnation  line  plots  of  N  and 
NO,  is  shown  next  to  the  wall,  in  refer 
ence  to  “no  plasma”  case.  Significant 
difference  is  noted  in  the  mass  frac 
tions  of  N  and  NO,  due  to  plasma 
DBD  actuator  in  plasmaS  configura 
tion.  For  N,  the  effect  is  seen  up  till 
25  mm  distance  away  from  wall. 


FIG.  25.  Effect  of  plasma  actuator  on 
N2  near  the  wall,  (a)  Effect  of  plasmaS 
configuration  on  stagnation  line  plot  of 
N2  near  the  wall,  and  its  (b)  zoom  in 
section.  A  significant  difference  is  seen 
in  curves  for  N2  up  to  30  mm  distance 
away  from  wall. 


recombining  into  N2,  O2,  and  NO.  For  the  distance  of  2  mm 
away  from  the  wall,  even  NO  species  vanishes  due  to  the 
super-catalytic  wall  boundary  condition.  Hence,  a  small  peak 
for  NO  is  observed  in  Fig.  24(b). 

Plot  of  N2,  in  Fig.  25,  shows  that  its  mass  fraction 
increases  in  two  phases  when  approaching  wall.  From 
30  mm  to  2  mm  distance  from  the  wall,  N2  monotonically 
increases,  after  which  there  is  a  dip  in  its  mass  fraction.  This 
is  due  to  mass  fraction  of  NO  increasing  at  a  distance  of 
2  mm  close  to  wall.  When  NO  start  decreasing  at  approxi¬ 
mately  1.5  mm  distance  away  from  wall,  N2  again  begins  to 
increase  in  its  mass  fraction.  Thus,  significant  effect  of 
plasma  DBD  actuators  is  noted  on  hypersonic  flows,  both 
with  and  without  thermo-chemistry. 

VIII.  CONCLUSION 

We  simulated  Mach  17  hypersonic  flow  over  a  cylinder, 
with  and  without  thermo-chemistry,  using  discontinuous 
Galerkin  method.  The  effect  of  micro- second  pulsed,  sinu¬ 
soidal  DBD  plasma  actuators  on  the  hypersonic  flow  was 
investigated  with  three  actuator  configurations.  A  significant 
change  in  the  surface  heat  transfer  was  predicted.  A  specific 
configuration  with  opposing  plasma  injected  momenta 
resulted  in  maximum  reduction  of  7.68%  in  the  stagnation 
point  heat  flux.  Low  magnitude  of  flow  velocity  near  the 
stagnation  point  makes  actuation  effective.  The  wall  jet 
induced  by  the  plasma  actuator  influences  the  temperature 
contours  near  the  actuator  location.  This  pull  and  push  mech¬ 
anism  is  responsible  for  the  increment  and  decrement  of 


existing  local  heat  flux  upstream  and  downstream  of  plasma 
actuator,  respectively. 

In  the  two  different  configurations  with  a  pair  of  plasma 
actuators  placed  at  a  distance  of  50  mm  (away  from  each 
other)  on  the  cylinder  surface,  the  second  configuration  gives 
overall  reduction  in  the  surface  heating  coefficient.  Thus,  we 
demonstrate,  that  the  micro- second  plasma  DBD  actuators, 
though  conventionally  thought  of  as  being  ineffective  for 
high  speed  flows,  do  have  significant  effect  on  the  surface 
heating  distribution  of  the  cylinder. 

In  addition,  the  effect  of  thermo-chemistry  on  the 
plasma  actuation  of  hypersonic  flow  was  also  studied,  using 

5  species  finite  rate  air  chemistry  with  17  reactions.  In  addi¬ 
tion  to  affecting  temperature  profiles  next  to  the  wall, 
plasma  actuators  also  affect  the  species  concentration  distri¬ 
bution  near  the  wall.  The  effect  is  significant  for  N2,  N,  and 
NO.  Based  on  the  effect  on  temperature  gradient  and  spe¬ 
cies  concentration,  a  maximum  reduction  of  4.04%  in  the 
stagnation  point  heat  flux  was  predicted  for  this  case.  This 
shows  promise  for  the  use  of  the  plasma  DBD  actuators  for 
the  integrated  surface  heating  reduction  in  hypersonic 
flows. 
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